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1. Inroduction 

The aim of this note is to prove a new discrepancy principle. The advantage of the new 
discrepancy principle compared with the known one consists of solving a minimization 
problem approximately, rather than exactly, and in the proof of a stability result. To 
explain this in more detail, let us recall the usual discrepancy principle, which can be 
stated as follows. Consider an operator eqution 

> . Au = f, (1) 

i-H ..... . 

Q\ . where A : H — > H is a bounded linear operator on a Hilbert space H, and assume that 
the range R(A) is not closed, so that problem (1) is ill-posed. Assume that / = Ay 
where y is the minimal-norm solution to (1), and that noisy data fg are given, such 
that ll/^ — JII < 5. One wants to construct a stable approximation to y, given fg. 
The variational regularization method for solving this problem consists of solving the 
^3 ■ minimization problem 

a 

F(u) := \ \Au - f s \ | 2 + e| \u\ \ 2 = min . (2) 



It is well known that problem (2) has a solution and this solution is unique (see e.g. [1]). 
Let ug >e solve (2). Consider the equation for finding e = e(5): 

\\Au s , e -f s \\ = C6, (3) 

where C = const > 1. Equation (3) is the usual discrepancy principle. One can prove 
that equation (3) determines e = e(5) uniquely, e(S) — > as S — > 0, and ug := us jt {S) V 
as 5 —>■ 0. This justifies the usual discrepancy principle for choosing the regularization 
parameter (see [1] and [2] for various justifications of this principle, [3] for the dynamical 
systems method for stable solution of equation (1), and [4] for a method of solving 
nonlinear ill-posed problems). 
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The drawback of this principle consists of the necessity to solve problem (2) exactly. 
The other drawback is the lack of information concerning stability of the solution to (3): 
if one solves (2) approximately in some sense, will the element us, e (S) ( with e(8) being an 
approximate solution to (3)) converge to y? 

Our aim is to formulate and justify a new discrepancy principle which deals with the 
both issues mentioned above. 

Our basic result is: 

Theorem 1. Assume: 

i) A is a bounded linear operator in a Hilbert space H , 

ii) equation Au = f is solvable, and y is its minimal-norm solution, 
iti) 11/(5 ~~ /II < \\fs\ \ > C5, where C > 1 is a constant. 

Then: 

j) equation ( 3) is solvable for e for any fixed 8 > 0, where us y€ is any element satisfying 
inequality F(us, e ) < m + (C 2 — 1 — b)5 2 , F(u) := \\A(u) — fs\\ 2 + e||-u|| 2 , m = m(8, e) := 
inf u F(u), b = const > 0, and C 2 > 1 + b, 

and 

jj) if e — e (^) solves (3), and u$ :— us t6 (5), then lim^o — y\ \ = 0. 
In Section 2 proof of Theorem 1 is given. 
2. Proof 

Proof of Theorem 1. Let us first prove the existence of a solution to (3). We 
claim that the function h(S,e) := HAu^ — f§\\ is greater than C5 for sufficiently large 
e, and smaller than C5 for sufficiently small e. If this is proved, then the continuity of 
h(8, e) with respect to e on (0, oo) implies that the equation h(5, e) = C5 has a solution. 

Let us prove the claim. As e — > oo, we use the inequality: 

e|K, e || 2 < F{u 5 ,e) < m + (C 2 - 1 - b)5 2 < F(0) + (C 2 - 1 - b)5 2 , 
and, as e — > 0, we use another inequality: 

||^-/,|| 2 < F(u 5 , £ ) < m +{C 2 -l-b)5 2 < F(y) + (C 2 -l-b)5 2 = e \\y\\ 2 + (C 2 -b)6 2 . 
This inequality implies 

h 2 {5,t) <t\\y\\ 2 + {C 2 -b)5 2 . 

As e — > oo, one gets \\us, e \\ < ^ — * 0, where c > is a constant depending on S. Thus, 
by the continuity of A, one obtains 

hm h{6,e) = \\A(p)-f s \\ = \\f s \\>C6. 

As e — > 0, one gets 

liminf/i(5,e) = liminf(e| \y\ \ 2 + (C 2 - b)5 2 ) 1/2 <C5. 

Therefore equation h(5, e) = C5 has a solution e = e(8) > 0. 
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Let us now prove that if u$ :— us t e(S), then 



lim | \us — y\ \ = 0. (4) 

From the estimate 

F(u 5 = | \Au 5 - f 5 \ | 2 + e\ \u s \\ 2 < C 2 5 2 + e\ \y\ | 2 , 
and from (3), it follows that 

IKII < IMI, ( 5 ) 

where w is any solution to (1). We will use (5) with w = y and w — U, where U is a 
solution to (1) constructed below, and y is a minimal- norm solution to (1). 

Thus, one may assume that ug — 1 U, and from (3) it follows that Au$ — > / as S — > 0. 

This implies, as we prove below (see (8)), that 

AC/ = /. (6) 
We also prove below that from (5) it follows that 

Km ||u, -tf ||=0, \\U\\<\\y\\. (7) 

o— >U 

The minimal norm solution to equation (1) is unique. Consequently, (7) implies U = y. 
Thus, (4) holds. 

Let us now prove (6). We have 

(/, v) = \im(Au 5 , v) = \im(u 5 , A*v) = (U, A*v) = (AU, v) W G H. (8) 

Since v is arbitrary, (8) implies (6). 

Finally, we prove (7). We have us — 1 U. Thus, \\U\\ < liminf^o ll^ll- Inequality (5) 
implies limsup^Q \\us\\ < \\U\\. Consequently, \\us\\ — > \\U\\. It is well known that the 
weak convergence together with convergence of the norms imply in a Hilbert space strong 
convergence. Therefore lim^o \ \us ~ U\ \ = 0. Taking w = y in (5), and then passing to 
the limit 5 — > in (5), yields inequality (7). Thus, both parts of (7) are established. 

Since U solves equation (1), and \ \U\ \ < \ \y\\, it follows that U = y, and (4) holds. 

Theorem 1 is proved □. 
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